Electronic structure and weak electron-phonon coupling in TaB2 
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We present electronic structure calculations together with resistivity, susceptibility, and specific 
heat measurements for TaB2 to search for the recently contradictorily reported superconductivity 
and to study related normal state properties. We ascribe the absence of superconductivity down to 
1.5 K for our TaE>2 samples to the generally weak electron phonon coupling derived from comparison 
of the calculated and measured specific heat constants. For the E2 3 and the Bi 3 V point phonons 
we derive from the calculated deformation potentials very small electron phonon couplings for these 
modes, opposite to the strong coupling of the E2 g mode in MgB2, probably responsible for its high 
T c . In comparison to MgB2, we discuss the origin of the quite different features in the density of 
states and of the Fermi surfaces. The differences are mainly due to the strong hybridization between 
Ta 5d and B 2p states outside the hexagonal basis plane. 

74.25-q,71.20.-b 



I. INTRODUCTION 



II. METHODS 



The recent discovery of superconductivity in MgB 2 0] 
has initiated an immediate broad research activity due to 
the surprisingly high transition temperature T c ~ 40 K in 
a seemingly ordinary s—p metal. Investigation of related 
s - p diborides MB 2 (M = Li, Be, Al) [§-[1 and a series 
of isostructural transition metal diborides (M = Sc, Ti, 
Zr, Hf, V, Ta, Cr, Mo, Nb) [§f) has shown that only few 
of them seem to be superconducting, and then only at 
very low temperatures. Since the understanding of the 
pairing mechanism in MgB2 is still in its early stages, 
a study of the isomorphic compounds with low transi- 
tion temperatures or with absence of superconductivity, 
respectively, could be helpful in clarifying the expected 
very substantial electronic differences. 

Contradictory reports about superconductivity in 
TaB 2 have appeared; one found superconductivity at 
T c ~ 9.5 K, |^J^] while another found no superconduc- 
tivity down to 4.4 K. @ TaB 2 is isovalent with NbB 2 , 
where the occurrence of superconductivity, or at least the 
value of T c , is likewise uncertain. NbB 2 was reported to 
be superconducting at 3.87 K, |j| superconducting only 
at 0.62 K, Q and not superconducting above 0.37 K. 
JTof Clearly there is sample dependence for both TaB 2 
and NbB 2 , and no doubt for many other diborides as 
well. 

In an attempt to begin to settle some of these discrep- 
ancies, we address TaB 2 specifically in this paper with a 
joint theoretical (Sees. II A and III A) and experimental 
(Sees. II B, III B) investigation. In Sec. Ill we provide an 
analysis of its electronic structure, and contrast it with 
that of MgB 2 . We also present specific heat, susceptibil- 
ity as well as resistivity data for two almost single phase 
samples. One of them exhibited a significant amount of 
boron vacancies. 



A. Band structure calculations 

We calculated the electronic structure of TaB 2 in the 
hexagonal space group (SG) PG^/mmc (No. 191) with the 
lattice constants a = 3.082 A and c = 3.243 A. ||| The 
frozen phonon calculations for the E 2s and Bi g modes 
where done in the orthorhombic SG Cmmm (No. 65) 
and in the trigonal SG P3ml (No. 164), respectively. 

Our band structure calculations were performed using 
the full-potential nonorthogonal local-orbital minimum- 
basis scheme (FPLO) Jll|] within the local density ap- 
proximation (LDA). In these scalar relativistic calcula- 
tions we used the exchange and correlation potential of 
Perdew and Zunger. [jl2j Ta 5s, 5p,6s, 6p, 5d states and 
B 2s, 2p, 3d, were chosen as minimum basis set for the 
valence states. All lower lying states were treated as core 
states. The inclusion of the relatively extended Ta 5s, 5p 
semi-core states as band states was done because of the 
considerable overlap of these states on nearest neighbors. 
This overlap would be neglected if they were treated as 
core states in our FPLO scheme. Accounting for this 
overlap is of importance especially for the calculations 
of phonon frequencies that we report. B 3d states were 
added to allow for boron polarizability. The spatial ex- 
tension of the basis orbitals, controlled by a confining 
potential [p"3| (r/ro) 4 , was optimized to minimize the to- 
tal energy. The self-consistent potentials were carried out 
on a /c-mesh of 20 fc-points in each direction of the Bril- 
louin zone, which corresponds to 481, 784 and 1221 in 
irreducible part for V63/ mmc, P3ml, and Cmmm SG, 
respectively. Our FPLO approach has been proved by 
recent comparison to FLAPW results in many cases to 
have an absolute accuracy of lmHartree/atom for the 



total energy and a much higher accuracy for the relative 
total energy changes. 

B. Sample preparation and measurements 

Two different TaB2 samples, one with near stoichio- 
metric and one with boron deficient compositions, were 
prepared starting from the pure elements Ta (> 99.9%) 
and B (> 99.9%) by arc melting furnace under purified 
argon atmosphere. To ensure a better homogeneity the 
samples were turned several times. The phase content 
was checked by X-ray diffraction (XRD) using Co-Ka 
radiation. 

The samples contain the hexagonal AIB2 phase mainly, 
95% in sample No. 1 and 97% in sample No. 2. In both 
samples small amounts (< 5%) of a second phase were 
found which are B for sample No. 1 and Ta3B4 for sample 
No. 2. The compositions of the main phases measured by 
electron probe microanalysis in the WDX modus showed 
a composition of about TaB2.o3 (sample No. 1) and com- 
positions of TaBi.29 (sample No. 2). 

Furthermore, the XRD resulted in quite different lat- 
tice parameters for the two samples ranging from a = 
3.098±0.002 nm and c = 3.224±0.003 A of the B-deficient 
sample to a = 3.067 ± 0.002 A and c = 3.286 ±0.006 A 
of the B-rich sample. 

The specific heat of both TaB2 samples was measured 
in the temperature range between 2 and 16 K using a 
Quantum Design PPMS relaxation calorimeter. The ad- 
denda which were determined in a separate run were sub- 
tracted in order to obtain the specific heat data for TaB2. 

Magnetization measurements have been performed us- 
ing a Quantum Design-SQUID magnetometer in the tem- 
perature range down to 1.8 K. Resistivity measurements 
down to 1.5 K have been performed using the standard 
four point method. 



electrons of Ta with respect to Mg, the Fermi level has 
shifted from the bonding B states below the hybridiza- 
tion gap in MgB2 to the anti-bonding states above this 
gap in TaB2. From the very similar shape of the partial 
B and Ta DOS, a strong hybridization between B 2p and 
Ta 5d states is obvious, mentioned below in more detail. 
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FIG. 1. Total DOS and partial DOS of TaB 2 . The upper 
panel gives the total DOS, the middle and the lower panels 
show the contribution of Ta and B, respectively. The Fermi 
level is at zero energy. 

The calculated value of the density of states at the 
Fermi level N(sf) is slightly higher for TaB2 (N(ep) 
= 0.91 states/ (eVx cell)) than for MgB 2 (M^e) = °- 71 
states/eVxcell), in agreement with Refs. |l^,[l7|]). This 
corresponds to a bare specific heat coefficient 70 = 2.14 
mJ/(molexK 2 ) for TaB 2 . 



III. RESULTS AND DISCUSSION 
A. Theoretical results 

In Figure 1 we display the total as well as the atom 
decomposed density of states (DOS) of TaB 2 . The B 2p 
and Ta 5d states share almost equally in the occupied 
valence bands in the region -10 to -2 eV (the Fermi level 
E = Ef is taken as the zero of energy). Our DOS is in 
agreement with that of P.P. Singh Jl4| ], which we became 
aware right after the completion of our study. 

A striking difference in comparison to MgB2 is the 
dominating contribution of Ta 5d states to the DOS at 
Fermi level, which contribute about 70% of the total 
DOS; in MgB 2 the DOS at Fermi level is dominated by 
B-2p states (see Figure 2 for comparison |l5j ) . Although 
a rigid band picture is very limited in this case, a va- 
lence analysis shows that due to the 3 additional valence 
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FIG. 2. Total DOS and partial DOS of MgB 2 . The partial 
DOS are calculated by projection to the minimum orbital 
basis. The Fermi level is at zero energy. 

Figure 3 shows the band structure of TaB 2 along the 
symmetry lines of the hexagonal cell. As already men- 



tioned above, the bonding B a states, which lie in the 
region -10 eV to -2 eV and are highlighted in the mid- 
dle panel, are completely filled. These states, which are 
unoccupied along the T-A direction in MgB2 (compare 
to Figures 1 of Refs. [ ^6p^ ]), lie now between -5 eV and 
-2 eV for this symmetry line and show almost 3 eV dis- 
persion along the hexagonal axis (r-A) compared to 0.6 
eV in MgB2. The two-dimensional character of these 
states in MgB2 is obviously destroyed in TaB2. Further- 
more, we find around the A-L-H plane of the fc-space a 
strong hybridisation of those B a states with the Ta 5d xz 
and 5d yz states, indicated in the lower panel of Figure 3. 
Once more we like to emphasize the difference to MgB2, 
where these states are of nearly pure B a character. 
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FIG. 3. Band structure and band characters of TaB2 . 
The middle and the lower panels show the band characters 
of the B 2p x , y orbitals and the Ta 5d X z, y z orbitals, respec- 
tively. The line width is scaled with the orbital weights of the 
corresponding orbitals. 

The corresponding Fermi surface (FS) of TaB2 is 
shown in Figure ^|. For convenience of comparison with 
MgB 2 (compare to Fig. 3 in Ref. |l7|), we have chosen the 
A point as the center of the hexagonal prism. All three 



sheets of the FS of TaB2 are electron-like. Because of 
the strong dispersion of the antibonding B a - Ta 5d XZjVZ 
states, they build closed FS around the A point (see mid- 
dle and lower panel of Fig. 4), where the hole- like quasi 
two-dimensional tubes are found in MgB2 . The large FS 
in the upper panel of Fig. [| is due to Ta 5d states, the 
contribution of B states to this sheet is almost negligible. 




FIG. 4. The three different sheets of the Fermi surface of 
TaB2, all sheets are electron-like. The A point corresponds to 
the center of the hexagonal prism, the F point is the midpoint 
of the lower and upper hexagon. 

To account for the experimental uncertainty in the lat- 
tice constants for different samples |18|] , we also investi- 
gated the influence of different lattice constants for the 
experimentally reported range on the electronic struc- 



ture. The changes for the relevant features in the band 
structure are negligible, the DOS is basically unchanged, 
N(ef) varies by less then 2%. 

For a rough estimate of the electron-phonon (el-ph) 
coupling in TaB2, we calculated the phonon frequen- 
cies and the deformation potential of the E2 9 (in-plane 
displacement of the borons) and the Bi ff (borons dis- 
placed along z in different directions) zone-center phonon 
modes. Their frequencies are 98 meV and 85 meV, re- 
spectively. For the corresponding frequencies in MgB2 
Kortus et al. JlTj reported 58 meV and 86 meV, respec- 
tively. For AIB2 118 meV and 60 meV, respectively, were 
calculated. |l^] Already from the strong hardening of the 
calculated E2 9 frequency compared with MgB2 one can 
conclude a strongly reduced electron phonon coupling of 
this mode. 

Figure 5 shows the calculated band structure for the 
frozen E2 ff phonon mode of TaB2 with a B displacement 
of Aub = 0.018 A. The B bond stretching mode splits 
the antibonding B <r-Ta 5d xz _ yz bands along the T-A line. 
For an averaged split Ae^ / Aw b ~9 eV/A we find a defor- 
mation potential De 2 ~ 4.5 eV/A about 3 times smaller 
than -D_e 2g in MgB2. Jl6| Calculating larger elongations 
up to the actual rms we found a nearly linear dependence 
of the deformation potential De 2 on the elongation Au. 




FIG. 5. Band structure for frozen E2 9 mode, plotted along 
the same lines as in Fig. |^ to facilitate comparison. The 
two-fold degenerate bands in Fig. |3| crossing the Fermi level 
along the F-A direction, are split due to the symmetry reduc- 
ing phonon mode. 

Following Ref . []l6f , we can estimate the coupling from 
this mode alone using Eq. (2.34) of Kahn and Allen 
fjfjf for the EP matrix element in terms of De 2q - With 
Ne 2b (£f) — 0.28, resulting from a summation of the cal- 
culated orbital projected DOS for the orbitals mainly in- 
volved in the deformed band, we find 



X E2g = N E2g (e F )[ 



0.055. (1) 



Due to the high frequency oj and the smaller deformation 
potential De 29 , we get a coupling smaller by a factor of 
about 20 compared to MgB2. 

From a corresponding calculation of the coupling con- 
stant for the Bi g mode, reported jl^] to be softened in 
AIB2 (lu = 60 meV) , we find an even slightly smaller con- 
tribution compared with the already weak coupled E2 ff 
phonon. 

Assuming similar coupling to the acoustic phonons as 
reported jllj] for MgB2 or AIB2 and comparable contri- 
butions of other modes, the total el-ph coupling constant 
might be no more than A ~ 0.2. 



B. Experimental results 

The results of the specific heat measurements are 
shown in Figure ^. Since no superconductivity was ob- 
served at least in the temperature range down to 1.8 K 
from magnetization measurements (see Figure |?]), and 
down to 1.5 K from resistance measurements (see Fig- 
ure |^) , the results of the specific heat measurements are 
shown in the c p /T vs. T 2 plot. Thus, assuming stan- 
dard normal metal behaviour, the Sommerfeld coefficient 
7 was determined using the relation 



12 



c p /T = 1 +^RTr i Q- 3 T 2 , 



(2) 



where R being the ideal gas constant and O the initial 
Debye temperature. 




The sum on j runs over the two B atoms of the E2 S 
mode, for Mb we used the isotope averaged mass of B. 
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FIG. 6. Specific heat c p (T)/T versus T 2 of two TaB 2 
samples at zero magnetic field. The lines correspond to linear 
regression fits of the experimental data for temperatures T > 
8 K (sample 1) and T > 4 K (sample 2), respectively. 

The investigated TaB2 samples show distinct anoma- 
lies at temperatures around 4 K which are, however, not 
due to superconductivity as was found from magnetiza- 
tion and resistance measurements. 

Sizable deviations from the linear fits, shown in Figure 
^ are observed especially for sample No. 1. Similar devi- 
ations as found for these two samples are known as well 
from MgB2 |2l| |23]|. Since such anomalous contributions 
do not essentially change in the superconducting state of 



Mgl$2 , we ascribe them to lattice effects related to boron 
disorder present in all diborides. 
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FIG. 7. Temperature dependent magnetization M of two 
TaB2 samples in the temperature range 1.8 K < T < 10 K at 
an applied magnetic field of H = 10 Oe. 
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FIG. 8. Temperature dependence of the resistivity nor- 
malized to its values P296X = 350 /if2cm (sample 1) and p29ai< 
= 199 /iflcm (sample 2) at T — 296 K in the temperature 
range 1.5 K < T < 296 K. 

Naturally, the initial Debye temperature ©o 417 ± 
2 K is smaller than the corresponding value reported for 
MgB 2 of 750 to 800 K H,||j2l). But it is somewhat 
harder than the value one might expect from a simple 
scaling with the square root of the total mass ratio: 369 
± 11 K. This is in line with the hardening calculated for 
the optical phonons at the 7-point as reported above. 

TABLE I. Experimental results: Specific heat (Sommer- 
feld constant 7, (column 2), Debye energy Oo (column 3), 
residual resistivity ratio RRR (column 4), and resistivity 
P296K at T = 296 K. 
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The theoretically estimated value of the el-ph coupling 
constant A is in accordance with those small A values de- 
rived by comparing the experimental value of the Som- 
merfeld constant (see figure ||) with the calculated DOS 
at the Fermi level N{ep) 



le Xp = (tt 2 /3)4(i + \)N(e F ) = 7o (l + A). 



(3) 



Using our calculated 70 = 2.14 mJ/molexK 2 and the 
measured j exp = 2.8 mJ/molexK 2 of sample No. 1, we 
obtain an empirical value of Xno. 1 = 0.3. If this esti- 
mate is correct, a sizable contribution of relatively low- 
frequency phonons involving Ta-vibrations can be ex- 
pected. 

The low value of j exp = 2 mJ/molexK 2 for sam- 
ple No. 2 can be ascribed to the significant number of 
boron vacancies in understoichiometric TaB2-a samples. 
Within a very crude estimate, supposing a rigid band 
behaviour and using the atomic partial density of states 
Nto,( e f) and Nb(sf) (see Figure 1) we assume for the 
effective bare specific heat constant 7o,e// 

7o,e// < (Tr 2 /3)k 2 B {N Ta {e F ) + (2 - 5)N B (e F )). (4) 



For TaBi.29 (sample No. 2) this yields 70, e// = 1- 
corresponding to Xno. 2 = 0.053. Considering addition- 
ally broadening effects in the DOS due to disorder in the 
vacancy distribution and other impurities, that A value 
can be taken as a lower bound. 

Adopting a standard value of the Coulomb pseudopo- 
tential fi* = 0.13 one arrives at negligible values of the 
transition temperature T c ~ 10~ 7 K, irrespectively of the 
details of the shape of the Eliashberg function, i.e. the av- 
eraged phonon frequency. If even the pseudo potential fj,* 
would be ignored, T c would not exceed 18 mK. Thus, the 
experimental evidence for the absence of superconductiv- 
ity, possibly even down to several 100 mK, || becomes 
very plausible already in the traditional electron-phonon 
Migdal-Eliashberg picture. 



IV. CONCLUSIONS 

Although TaB 2 occurs in the same crystal structure 
as MgB2, it should not be considered as a close relative 
of MgB2 with respect to the electronic structure (even 
if it would be found superconducting at very low tem- 
peratures). The major differences occur due to: (i) the 
different band filling because of the three additional va- 
lence electrons of Ta with respect to Mg, resulting in a 
shift of the the Fermi level from the bonding B-cr states 
in MgB2 to the antibonding B-cr - Ta-<i hybrid states in 
TaB2, (ii) the strong out of plane hybridization of the 
B 2p states with Ta, (iii) the weak electron phonon cou- 
pling, especially of the E2 9 mode, which is strongly cou- 
pled in MgB2- 

According to our experimental results, TaB2 is not su- 
perconducting down to 1.5 K. Thus, we confirm the ear- 
lier results of Refs. pQ,p| and disprove at the same time 



the speculation about superconductivity around 9.5 K 
reported in Ref. Because in our opinion the main 
reason for the absence of superconductivity in TaE?2 is 
the different position of the Fermi level (with respect to 
MgB2), huge hole doping might "reintroduce" supercon- 
ductivity at relatively high temperatures. 

In other words, the results obtained here suggest that 
the empirical absence or low-temperature superconduc- 
tivity established in many transition metal (or rare earth) 
diborides with electrons as t he p otentially paired charge 
carriers, stressed by Hirsch |2q| , might be explained in 
the traditional electron-phonon picture simply by a weak 
electron-phonon interaction. In that case there is no 
need to explain this behavior by the absence of a sophis- 
ticated Coulomb interaction driven non-phonon mecha- 
nism which works exclusively for holes [ p5[ . 

Corresponding studies for other related transition 
metal diborides of experimental interest will be published 
elsewhere. 
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